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ABSTRACT 

Context. The existence of time-energy correlations in flare occurrence is still an open and much debated problem. 

Aims. This study addresses the question whether statistically significant correlations are present between energies of successive flares 

as well as energies and waiting times. 

Methods. We analyze the GOES catalog with a statistical approach based on the comparison of the real catalog with a reshuffled one 
where energies are decorrelated. This analysis reduces the effect of background activity and is able to reveal the role of obscuration. 
Results. We show the existence of non-trivial correlations between waiting times and energies, as well as between energies of 
subsequent flares. More precisely, we find that flares close in time tend to have the second event with large energy. Moreover, after 
large flares the flaring rate significantly increases, together with the probability of other large flares. 

Conclusions. Results suggest that correlations between energies and waiting times are a physical property and not an effect of 
obscuration. These findings could give important information on the mechanisms for energy storage and release in the solar corona. 
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1. Introduction 

■ Solar flares are violent explosions of magnetic energy in the 
solar corona. The theory of magnetic reconnection (see Priest 
& Forbes 120021 for a review) represents the most plausible and 

I widely accepted explanation for flare occuiTence, although a def- 
inite and clear explanation of the mechanisms at the basis of flare 
triggering is still lacking. A well-established property of flare 
occurrence is the power law decay of the peak-flux energy dis- 

■ tribut ion (Lee et al |1993l Aschwanden et al l 19981 Crosby et al 
' I1998I I. a property shared by other stochastic physical phenom- 
ena like earthquakes (de Arcangelis et al 2006). Several mech- 
anisms have been proposed to reproduce the above experimen- 
tal findings. In the Rosner and Vaiana (RV) model (Rosner & 
Vaiana 1978 ) it is assumed that flare occurrence is an uncor- 
rected Poisson process where the energy grows at a rate pro- 
portional to the internal energy of the system. Given this as- 

i sumption, the energy storage follows an exponential temporal 
growth, interrupted at random when the system releases all the 
stored energy in a flare. This mechanism accounts for the power 
law in the size distribution and predicts correlations between the 
storage time and the released energy. More precisely, the later a 
flare will occur the higher the energy will be. Also the avalanche 
model (Lu & Hamilton [19911 Hamon et al |2002] |. which de- 
scribes solar flares as energy relaxation events in a system driven 
at a constant rate, coiTectly predicts scale-free behavior for the 
peak-flux energy distribution. Within this approach, flare occur- 
rence is again a Poisson process, as in the RV model, but occur- 
rence times and energies are uncorrected like in standard self- 
organized models (Jensen 1998 ). Finally, power law behavior for 
the size distribution is also consistent with models assuming cor- 
relations between the waiting times At of successive bursts, as in 
a shell model of magneto-hydrodynamic turbulence (BofFetta et 
al ll999l l. This model reproduces the experimental power law de- 



cay of the waiting time distribution P(At) which can also be ob- 
tained however by means of an uncorrected piecewise Poisson 
process (Wheatland et al l2002l Norman et al l200TT l. 

The investigation of correlations between flare energies and 
waiting times is a useful tool to distinguish among different trig- 
gering mechanisms. This problem has been addressed in a series 
of papers (Crosby et al 1998. Wheatland et al 1998 , Wheatland 
2000), providing no clear evidence for time-energy correlations. 
In particular, Wheatland et al. (1 19981 ) detected small correlations 
at short time and interpreted them as a spurious effect of obscu- 
ration, i.e. long-lasting events may hide subsequent small events. 
By properly taking into account obscuration effects, Wheatland 
concluded that no significant correlations are present. In this pa- 
per we present a novel analysis of experimental catalogs provid- 
ing evidence for time-energy correlations in flare occurrence. 



2. Data analysis and methods 

In our analysis we used data from the GOES 
(Geostationary Operational Environmental Satellite) catalog 
(ftp://ftp.ngdc.noaa.gov/STP/SOLARJDATA/SOLAR_FLARES) 
for the solar cycles 21, 22 and 23 in the years 1977 to 2007. We 
include only CI. 4 class flares in the analysis (peak flux energy 
larger than £0 = 1.4 x 10"^ Wm"^). Each flare in the catalog is 
characterized by its starting time f, and its peak-flux energy 
in the following called "energy" to simplify notation and mea- 
sured in unit of Wm"^, counting 38575 flares with E > Eq. In 
order to enhance the eventual presence of temporal and energy 
correlations among flares, we performed a recently proposed 
(Lippiello et al 2008a) statistical analysis of the catalog. We 
computed the probability to have a given quantity x, greater than 
X conditioned to Af, being smaller than T, P(xi > X\Ati < T). 
Here, Af, = f, - r,_i is the waiting time and the quantity x,- rep- 
resents, depending on cases, the energy of the i-th flare x, - Ei, 
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Fig. 1. Disti-ibution of P*{Ei > E\Ati < T) for E = 3Eo, 
T - Ih and 10^ realizations of the reshuffled catalog repre- 
sented by black circles. The continuous red line is the fit with 
a Gaussian distribution with an average value of 0.3619 (dotted 
line) and a standard deviation o-{E, T) - 6.5 10"^. The dashed 
line indicates the value P(Ei > E\Ati < T) - 0.4098 obtained in 
the real catalog. 



the energy of the previous flare jc, = £■,_! or the the energy 
ratio between successive flares, x, - Ei/Ei^i. The quantity 
P{xi > X\Ati < T) is given by P{xi > X|Af,- <T) = where 
N{X, T) is the number of couples of subsequent events fulfilling 
both X, > X and Af,- < T, whereas N{T) is the number of couples 
satisfying only the condition Af,- < T . In order to detect eventual 
correlations in the catalog, we introduced the probability 
diff'erence 6P(xi > X|Af, <T) = P{xi > X|Af, < T) - P(xi > X), 
where f (jc, > X) is the unconditional probability to have x/ > X. 
A positive (negative) 6P indicates a higher (lower) probability 
to find flares with x, > X if one restricts the analysis to couples 
with Af, > T. Hence, values of SP diff'erent from zero suggest 
the existence of correlations between x,- and Af,. Evidently, 
because of statistical fluctuations, 6P is never exactly equal to 
zero even in catalogs where xt and Af,- are uncorrected. In order 
to explicitly take into account the role of statistical fluctuations, 
we computed the quantity P*(x, > X\Atj < T), defined as P, 
but in a catalog where flare energies were randomly reshuffled. 
The comparison with thereshuffled catalog is the basis of the 
surrogate data technique (Theiller et al. 119921 Schreiber & 
Schmitz 2000), recently applied to investigate coherent struc- 
tures in space plasma (Sahraoui 1). In the reshuffled catalog x, 
was by construction uncorrelated to Af, and P*ixi > X\Ati < T) 
fluctuated around its average value P(x,- > X). The amplitude of 
these fluctuations defined the significance level, a-{X, T), which 
allowed us to distiguish between the presence and absence of 
correlations. The method is schematically presented in Fig.l for 
jc, =£■,., X = £■ = 3Eo and r = 1 h. f *(£, > ^lAf, < T) takes 
different values for each realization of the reshuffled catalog. 
We produced 10^ independent realizations of the catalog with 
reshuffled energies and observed that P*(Ei > E\Ati < T) is 
Gaussian-distributed with a mean P{Ei > E) and a standard 
deviation cr{E, T). Similar results were obtained for other 
values of E and T and for the other definitions of x, . Therefore, 
\6P(xi > X\Ati < r)| > cr(X, T) indicates that energies in the real 
catalog follow a significantly different organization than events 
in the reshuffled catalog. 




Fig. 2. ProbabiHty diff^erence 6P(Ei > E\Ati < T) (a) and 
5P(Ei^i > E\Ati < T) (b) as a function of E/Eq for different 
values of T. For each data point the error bar is the standard 
deviation o-(E, T). 




Fig. 3. (a) Probability difference 6P(Ei-i > E\Ati < T) as a func- 
tion of £■/£■() for r = Ih and different k values, (b) Same quantity 
as function of Ele,h for T - \\\ and different values of e,i,. For 
each data point the error bar is the standard deviation o-(E, T). 



3. Results 

In Fig.2a we plot 5P{Ei > E\Ati < T) and in Fig.2b 5P(£',_i > 
E\Ati < T) for different values of T. The standard deviation for 
each data point cr{E, T) is represented as the error bar In Fig. 2a 
we notice that SP{Ei > E\Ati < T) takes always positive values 
beyond error bars. This implies that in the real catalog the num- 
ber of couples fulfilling both conditions is greater than in cata- 
logs where energies and intertimes are uncorrelated. More pre- 
cisely we find that for each given value of E, SP{Ei > E\Ati < T) 
decreases by increasing T. This implies that the probability to 
find flare couples with the second flare energy higher than E de- 
creases if one includes events with larger Af in the analysis. This 
result is disagrees with what is expected according to the RV 
model, which predicts larger flares after longer waiting times. 
On the other hand, obscuration can represent a possible explana- 
tion of the above result. Indeed, according to the selection pro- 
cedure, for flares close in time, the second event is recorded in 
the catalog only if it produces an increase in the flux of at least 
40% of the level of the previous flare. Moreover, the light curve 
of large flares may hide smaller flares close in time. Obscuration 
effects then reduce the probability to find a small flare after a 
short waiting time and may introduce spurious correlations be- 
tween waiting time and the successive flare energy. 
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Fig. 4. Probability difference ^fC^,-! > E\Ati < T) as a function 
of E/e,i, for T - Ih (a) and T - IQh (b) for different values of 
en,. For each data point the error bar is the standard deviation 
cr(E, T). 




Fig. 5. Probability density difference SpiEi^i - E\ai < Af, < 
ai) versus EjE^ for different values of a\ and q'2- 



We then explicitly explored the role of obscuration on 
6P(Ei > E\Atj < T). As a first analysis we considered only flare 
couples with a temporal distance greater than k times the dura- 
tion of the first flare. Spurious effects due to obscuration should 
disappear by increasing k. We observed (Fig. 3a) that this extra 
condition did not reduce correlations between waiting times and 
energies, since they are still present also for the highest values 
of k. To further prove that correlations cannot be attributed to 
spurious effects related to obscuration, we performed the same 
analysis as in Fig. la, setting different lower energy thresholds 
e,i,. Indeed, by increasing the value of the minimum energy re- 
quired for flares to be included in the analysis, the percentage of 
flares hidden by obscuration was reduced and obscuration effects 
should progressively decrease. In Fig. 3b we plot the results for 
different energy thresholds e,/, as a function of E/e,/,. No signif- 
icant dependence on e,/, is detected, indicating that coiTelations 
cannot be attributed to obscuration effects. We conclude that for 
couples of events close in time it is more probable to have the 
second event with a high energy than for events distant in time. 

In Fig. 2b we performed the same analysis as in Fig. 2a for 
the correlations between the waiting time and the previous flare 
energy. In contrast to Fig. 2a the quantity 6P(Ei-i > E\Ati < T) 
is not always positive now. In particular for T = Ih, 5P(£'/_i > 



E\Ati < T) is positive at small E and negative at large E. In 
order to verify if the behavior at small T is affected by obscu- 
ration, we performed the same analysis as in Fig. 3b, imposing 
the condition on the lower energy threshold e,/,. In Fig. 4a we 
present 5P(£',_i > E\Ati < T) for T = Ih and different energy 
thresholds. We found that data strongly depend on e,/,, in partic- 
ular for e,/, = lOEi), 6P{Ei-\ > E\Ati < T) is comparable with 
the statistical fluctuations. This analysis suggests that the corre- 
lations observed at small At are related to obscuration effects. 
Conversely, data for T = 1 Oh are not significantly affected by 
e,/, and 6P{Ei-i > E\Atj < T) always takes positive values be- 
yond eiTor bars (Fig. 4b). This result indicates that coiTelations 
detected for events with temporal distances between Ih and lOh 
are not significantly affected by obscuration. We conclude that 
couples of events distant in time between Ih and lOh tend to 
have the first event with a higher energy than the one expected 
for a process where time and energy are uncorellated. 

To obtain more detailed information on coiTelations be- 
tween £■,_! and Af, we considered the quantity dp{E,a\,a2) = 
5P{Ei-\ - E\ai < Af; < ai), i.e. the difference between condi- 
tional and unconditional probability densities to have £■,_! = E 
if Af; e]Qfi,Q'2[- This quantity can be obtained deriving the cu- 
mulated probability 6P{Ei^\ > E\ai < Af,- < Q'2) with respect 
to E. Being a derivative, it is affected by higher fluctuations. For 
every choice of ai, the value of a2 was chosen to have the same 
number of events, about 9000, for all curves. We noticed that for 
A/ > ai = 20h the difference in the probability density dp was 
very close to zero, indicating the absence of correlations. 

We now discuss the results of Fig. 5 at fixed E/Eq and for 
varying ai. We first observed that for small of E, dp is always 
negative for all values of ai. The negative region is restricted 
to very small values of £ < I.SEq for ai - and increases to 
E < 23Eq for larger a\ . Within this interval dp is non monotonic 
in Q-i , in particular, it reaches the minimum value for o-i = 2h and 
then tends to zero. For Af, in a given range, negative values of 
dp imply that it is less probable to have a flare with an energy 
E, - E compared to the case where £, and Af, are uncorrelated. 
Results in Fig. 5 then imply that it is improbable to have small 
Af after very small flares and that these anticorrelations vanish 
for large Af. In particular the highest negative values of dp are 
obtained for At e {2h, Ah} after an £ < 2Eq flare. 

If these anticorrelations were a spurious effect of obscura- 
tion, we should rather observe a different behavior. Indeed, en- 
ergy reshuffling produces more events shortly after a large flare 
compared to the real catalog. Consequently, we should observe 
an excess, and not a deficit, of events after a small flare in the 
real catalog compared to the reshuffled one because the num- 
ber of events is conserved. We now discuss higher values of 
E > \OEq. In this case dp moves from negative values for ai -Q 
and a\ - Ih to positive values for ai - 2h and a\ - 4h and then 
tends to zero for higher values of ai . This indicates that there is 
a deficit of events with Af < 2h after energetic flares, whereas 
for Af €[2h,4h] the number of observed flares is larger than the 
one expected for an uncorrelated process. The above results sug- 
gest that after a large flare a recovery time of about 2 hours takes 
place, when only a small number of flares is observed. After this 
time the number of observed flares reaches a maximum value for 
Af < 4h and then decreases to the background level for larger Af . 
According to the dependence on e,/, observed in Fig.4a, the re- 
covery time can be a spurious effect related to obscuration. 

The number of events n(f) occurring at the time f after an 
energetic flare can be explicitly evaluated following the method 
of Lippiello (2008b). More precisely, we define as a "main" flare 
an event with the energy E > Emain - lOZiQ. Indicating with 
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t(hours) 



Fig. 6. (a) Average number of flares n(f) as a function of the time 
t occurring after a main flare with an energy E > £„,„,„ - IQEq. 
The data are averaged over 4057 main flares, (b) The difference 
Sp(E) between the energy distribution for the first m - 20 flares 
occurring after an event with a higher energy than £, and the 
entire catalog energy distribution and different Er values. 



f, the occurrence times of main-flares we compute the quantity 
N(t, Emain) — Zi "(^ ~ ti)®{ti+\ - > whcrc ®{x) is the Heaviside 
step function, and the sum extends over all main-flares. Here 
n(f - ?,) is the number of flares with an energy lower than lOZio 
occurring during the time t - ti after the i-th main flare, and the 
sum extends over all main-flares in the catalog. Assuming that 
n{t - ti) is time translationally invariant we obtain 



n(f) = 



(1) 



which represents the average number of flares occurring in a pe- 
riod of a duration f after a mainflare. The data (Fig. 6a) agree with 
the previous results, with a maximum number of flares at f ^ 4h, 
decaying at longer times. Notice that the asymptotic decay is 
consistent with a power law f with p - I, reminiscent of the 
Omori law for seismic sequences (Omori [T894| de Arcangelis et 
al l2006l) . 

We now discuss the existence of correlations between ener- 
gies of subsequent flares. In order to do so, we computed the en- 
ergy distribution p{E) of the first m flares occurring after a flare 
with a larger energy than a reference value E,-. Since the analysis 
is restricted to all flares with E > Eq, obviously p{E) coincides 
with the flare energy distribution of the whole catalog Pt{E), for 
Er = Eq. In Fig. 6b we plot 6p(E) - p(E) - priE) for m = 20 
and different values of E^. Similar results are obtained for other 
values of m. Deviations from priE) become more and more ev- 
ident for increasing values of Er- In particular, the higher Er is, 
the larger is the probability to have subsequent flares with higher 
energy and the lower is the probability to have small flares. 

In order to get further insights in the correlations between 
flare energies, we computed the quantity 6p{Ei - AEi^ilAtj < 
T), i.e. the difference between the conditional and the uncon- 
ditional probability density to have the energy of the subsequent 
flare A times the previous one. We found (Fig. 7a) that for T = Ih, 
6p(Ei = AEi-i\Ati < T) is significantly different from zero for 
all A values. This indicates that the energies of two subsequent 
flares are correlated. Moreover we observed that these correla- 
tions depend on the time separation between the two flares and 
are practically zero for Af, > lOh (Fig.7a). We verified that this 




Fig. 7. Probability density difference dp{Ei - AEi-[\Ati < T) as 
a function of A for different T (a) and for T - Ih and different 
eth (b). 

result is not affected by obscuration performing the same analy- 
sis for r = Ih and different lower energy thresholds e,/, (Fig. 7b). 
We found that curves for different e,/, coincide within statistical 
fluctuations. This indicates that energy correlations and their de- 
pendence on time separation are a physical property, not a spu- 
rious effect due to obscuration. Curves present a maximum for 
A > I, indicating that it is more probable to find the next flare 
with an energy close to but slightly higher than the previous one. 

4. Conclusions 

In conclusion, we presented a statistical analysis of the GOES 
catalog indicating the existence of time-energy correlations be- 
tween successive events not to be attributed to obscuration ef- 
fects. More precisely, we observed that for couples of events 
close in time (T < Ih), the second event tends to have a high 
energy. Moreover couples of events distant in time between Ih 
and lOh have the first event with a high energy. The analysis of 
the rate decay after large flares shows evidence that the largest 
number of events is detected about 4h after the occurrence of the 
main event. Finally the distribution of flare energies confirms 
that the higher the flare energy, the larger the number of sub- 
sequent events with high energy. The existence of time-energy 
correlations suggests the possibility of scaling laws relating time 
with the energy released in a flare. This is a still open question, 
which could provide interesting insights in the energy storage 
and release mechanisms at the origin of solar flare occurrence. 
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